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Abstract 

We study the influence of the surface orientation on the wetting behavior of bcc binary 
alloys, using a semi-infinite lattice model equivalent to a nearest-neighbor Ising antiferro- 
magnet in an external magnetic field. This model describes alloys that exhibit a continuous 
B2-A2 order-disorder transition, such as FeAl or FeCo. For symmetry-breaking surfaces like 
(100) an effective ordering surface field gi 7^ emerges. Such a field does not only crucially 
affect the surface critical behavior at bulk criticality, but also gives rise to wetting transitions 
below the critical temperature T c . Starting from the mean-field theory for the lattice model 
and making a continuum approximation, a suitable Ginzburg-Landau model is derived. Ex- 
plicit results for the dependence of its parameters (e.g., of gi) on the microscopic interaction 
constants are obtained. Utilizing these in conjunction with Landau theory, the wetting phase 
diagram is calculated. 
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1 Introduction 



Surface critical behavior at bulk critical points can be divided into distinct universality classes 
0. For a given bulk universality class, only gross surface properties determine which surface 
universality class applies, such as: whether or not the surface interactions exceed or are equal to 
a certain critical enhancement, or whether a surface field g\ coupling to the local order parameter 
exists. Recently it has been shown that the universal critical behavior at the surface of a bcc 
Ising antifcrromagnet and of a binary alloy undergoing a continuous order-disorder bulk transition 
depends crucially on the orientation of the surface with respect to the crystal axes |). The 
basic mechanism underlying this intriguing behavior is the interplay between broken translational 
invariance perpendicular to the surface and the symmetry with respect to sublattice ordering. For 
certain "symmetry-breaking" orientations an "effective" ordering surface field gi ^ emerges, 
which depends on physical parameters like temperature and bulk composition of the alloy. That 
such a field exists has already been pointed out in j|] in order to explain the persistence of 
surface order at a (100) surface above the bulk critical temperature T c , detected in a Monte Carlo 
simulation for the B2-A2 order-disorder transition in Fe-Al. 

The situation encountered for symmetry-breaking surfaces closely resembles the critical ad- 
sorption of fluids, where generically g\ ^ ||. However, in that case the microscopic origin of gi 
is quite different: It is an external field reflecting, e. g., the preference of the wall for one of the 
two components of the binary liquid mixture. The transition that takes place at the surface of 
the system in the presence of a field g\ ^ on approaching the bulk critical point has been called 
normal in If <?i = (and the surface interactions are not too strongly enhanced), another 
transition, called ordinary, occurs. In accordance with the fact that gi is a relevant scaling field, 
the ordinary and normal transitions represent different surface universality classes. 

In Refs. J2j and || the focus has been on the behavior at T — T c and a clear identification 
of the normal transition, which may also be regarded as a critical point wetting phenomenon . 
However, since g\ generally stays nonzero away from T c for symmetry-breaking surfaces, a variety 
of wetting phenomena may occur for T < T c . Below we will determine the wetting phase diagram 
for a (100) surface within the mean-field approximation, utilizing the continuum model derived 
in H . Our work complements previous studies on wetting in fee Ising antiferromagnets or binary 
alloys Q as well as on interface roughening at an antiphase boundary in the [100] direction in 
bcc binary alloys Q. 

The organization of the paper is as follows. In the next section we define our model, explain 
the difference between symmetry-breaking and symmetry-preserving surfaces, and then briefly 
discuss the discrete mean-field equations (Sec. ||). In Sec. || we introduce the Ginzburg-Landau 
model for the (100) surface derived in ||. This is then used in Sec. [I] to determine the wetting 
phase diagram. 



2 Lattice model 
2.1 Definition 

To model the continuous B2-A2 order-disorder transition in the binary (AB) alloys FeAl or 
FeCo, we consider a bcc Ising antiferromagnet with nearest-neighbor (NN) interactions of strength 
J < . The spin variable o~i takes the values +1 or —1 depending on whether lattice site i is 
occupied by an A or B atom. Within the grand-canonical ensemble, the Ising Hamiltonian reads: 

i zesurf 

where Yl/ij) runs over all NN bonds. The bulk field H serves to adjust the composition of 
the alloy and represents the chemical potential difference between A and B atoms. A nonzero 
surface field Hi occurs generically in binary alloys, giving rise to surface segregation effects. One 
has Hi > (Hi < 0) if A (B) atoms tend to segregate at the surface. It is important not to 
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confuse Hi with an "ordering" field, which couples directly to the local order parameter. Only if 
Hi distinguishes one of the two sublattices at the surface, as is the case for symmetry-breaking 
surface orientations where the surface sites belong to a single sublattice only (see below), will Hi 
contribute to an "effective" ordering field gi . But even then Hi ^ is not a necessary condition 
for gi (cf. Sec. g). 

The average concentration (or occupation probability) Cj of A atoms on lattice site i can be 
written in terms of the mean magnetization m, = (o~i) of spin o~i as = ^(1 + m,j). In the 
ordered (52) phase, the bcc lattice splits into two interpenetrating sc sublattices a and f3 with 
bulk magnetizations ^ m^, which are preferentially occupied by A and B atoms, respectively 
(cf. Fig. 1). The disordered (A2) structure is characterized by = =: mdis- The bulk order 
parameter is defined by 

fa = \ (< ~ mfj . (2) 

2.2 Symmetry properties of the surface 

Let us define more precisely what is meant by symmetry-breaking and symmetry-preserving 
surfaces. Consider a uniform translation r a p of the crystal lattice that maps a-sites into /3-sites. 
In an infinite system without free surfaces or a finite system with periodic boundary conditions, 
the Hamiltonian TL = 7i{<x;}, Eq. (|l|), is invariant under the transformation: 

°i ^ a 'i = a i+T a p ■ (3) 

This symmetry is spontaneously broken below T c where the mean magnetizations m? and 
transform into each other under (g), so that 0b — > —(f>b- For a system with a free surface, the 
invariance of the Hamiltonian still holds if T a p can be taken parallel to the surface as is the case 
for the (110) orientation (Fig. 1). Then the surface is called symmetry-preserving. One may 
convince oneself that a surface with Miller indices (711712773) is symmetry-preserving if and only 
if rii + n 2 + n 3 is even. (By convention, we use the cubic unit cell of the bcc lattice here.) The 
order parameter, which becomes a local quantity <p — <p n depending on the discrete layer index n, 
vanishes identically above T c since the <f> — ► —(f) symmetry of the bulk system is retained and no 
enhanced surface interactions have been admitted in ([!]). By contrast, if either one of the lattice 
planes parallel to the surface belongs to a single sublattice, as is the case for the (100) surface (Fig. 
1), no translation T a p parallel to the surface exists and the surface is called symmetry-breaking. 
Generally, a (711712^13) surface is symmetry-breaking if n\ +712 + 773 is odd. Then the Hamiltonian 
is no longer invariant under (g) but changes by an amount proportional to the total number of 
surface sites. (One may again consider a finite system but impose periodic boundary conditions 
only in the directions parallel to the surface.) Thus the <j) — > —<j) symmetry of the bulk system is 
generically broken, and the order parameter will be nonvanishing at least locally near the surface 
even if T > T c . 

The symmetry properties of the surface must also show up in the context of suitable continuum 
(Ginzburg-Landau) models. Let 4> s be the value of the order parameter at the surface. For 
symmetry-preserving surfaces, the surface contribution to the Landau free energy will only contain 
even powers of (f> s . By contrast, arbitrary odd powers are expected to occur for symmetry-breaking 
surfaces due to the loss of the <f> — > — symmetry. In particular, the coefficient of the linear 
term may be identified with an "effective" ordering surface field gi =/= 0. Of course, in order to 
estimate the magnitude of 171 and its dependence on physical parameters such as temperature and 
bulk composition, the parameters of the continuum model must be related explicitly to lattice 
quantities (see Sec. ||). 

2.3 Mean-field (Bragg- Williams) approximation 

Owing to the spatial inhomogeneity along the z axis perpendicular to the surface, exact 
treatments of the model (|l|) are very hard, and one usually has to rely on approximate techniques 
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such as the mean-field (MF) or Bragg- Williams approximation. The MF equations read (with 
ks = Boltzmann's constant): 



= tanh 



(i) 



(4) 



where Hi = H , Hi = H + Hi for bulk and surface sites, respectively. The sum Yy runs over all 
NN sites of i. The mean magnetizations rrii vary with the index n = 1, 2, . . . labeling the lattice 
planes along the z axis, but are the same on each sublattice within a layer. Thus for the (110) 
surface, two variables are needed to describe the state of each layer: 



nii — to" for i € layer n, subl. a , TOj = to^ for i G layer n, subl. (3 . 
The local order parameter 4> n is conveniently defined as 
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to, 



to;. 



For the (100) surface, one may write 



for i £ layer n , 



(5) 



(G) 



(7) 



since each lattice plane belongs to a single sublattice. The definition of the local order parameter 



requires more care. The obvious choice 



;(— l) n+1 (m„ — rn n+ i) is physically reasonable but 



causes considerable problems in the continuum limit, as explained in detail in ||. It is favorable 
to adopt the more symmetric definition 



;(-!) 



n+l 



■(to„ + i + m„_i) 



(8) 



which treats the preceding and succeeding layer on an equal footing. 

The MF equations (||) have been studied in || both for (110) and (100) surfaces via the 
"nonlinear-mapping" technique [|oj. It has been shown that in the case of the (100) surface a 
nonvanishing order parameter profile appears for T > T c . The characteristic length scale that 
governs its exponential decay at large z diverges as T — * T c . Precisely at T — T c , the decay takes 
the form z~Pl v (where /3 = v = 1/2 within MF theory), which is one of the signatures of the 
normal transition. In the case of the (110) surface, the local order parameter vanishes for T > T c 
since to" = to^. Nevertheless one obtains a nontrivial magnetization profile, so that A rich and 
B rich layers alternate as one moves along the z direction. However, the length scale associated 
with this profile remains finite at T = T c . 



3 Continuum (Ginzburg-Landau) model 

In proceeding to a suitable continuum description one must be aware of several novel features 
arising for Ising antifcrromagnets or binary alloys, which are not present in simpler systems 
equivalent to the Ising ferromagnet. First, the Landau expansion of the surface free energy should 
look different for distinct orientations, due to the loss of the 4> — > ~4> symmetry for symmetry- 



breaking surfaces (cf. Sec. 2/2). Second, one has to take into account so-called "non-ordering" 
(or noncritical) densities. In the alloy picture these are needed to account for spatially varying 
profiles of, e. g., the local concentration that could not be described by the order parameter alone. 

Non-ordering densities introduce additional length scales that may compete with the order 
parameter correlation length £{,. It has been shown in a study of wetting in fee Ising antifer- 
romagnets that this competition may even lead to nonuniversal critical wetting exponents, and 
that the number of non-ordering densities depends on the orientation of the surface. Of course, 
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noncritical densities do not affect the asymptotic surface critical behavior if the bulk transition is 
continuous (as in our case of a bcc antiferromagnet), since then the diverging correlation length 
£b dominates all other length scales. However, we are interested here in phenomena below T c , 
where is finite. Accordingly, length scales associated with such spatially varying non-ordering 
densities may well be of the same order and important. 

To become more specific, let us recall the Ginzburg-Landau model derived and critically 
examined in Ref. || for the case of the symmetry-breaking (100) surface. This is based on a 
free-energy functional of the form 



dz 



+ hW)\\ + Ufa) , (9) 



where <f> s = 4>(0). The Landau expansions of the bulk and surface free-energy densities read 

M<t>) = \<P + + o{<f) , /.OM = -$i& + !# + 0($). (io) 

As expected, arbitrary odd powers of <f> s occur in the expansion of f s (<p s ). In view of the above 
discussion it is remarkable that no spatially varying nonordering density appears in (||). The 
reason is that lattice planes belonging to sublattice a and (3 alternate along the [100] (or z) 
direction. Hence the order parameter profile uniquely determines the segregation profile and 
vice versa. By contrast, the nontrivial segregation profile present for the (110) surface above 
T c provides a typical example of a spatially varying non-ordering density i/j appearing in the 



corresponding free-energy functional 11 . The "segregation field" Hi then couples linearly to ip, 
but no terms linear in the local order parameter and thus no ordering surface field are present. 

Deriving the continuum theory from the lattice model has the virtue of relating the "phe- 
nomenological" coefficients in ([l0]) explicitly to the "microscopic" parameters: 

k B T ' 4|J| ' 1 4|J| ' V ; 

which are dimensionless measures of the spin coupling strength and the bulk and surface magnetic 
fields. The "bulk" coefficients a, b, and c are independent of surface properties and depend on 
K and h only: a = a{K, h) etc. As usual, a varies linearly with the reduced temperature 
t = (T — T c )/T c near T c (at fixed magnetic field h), whereas b and c are positive constants to 
lowest order in t (see ||). The "surface" parameters are found to be: 

\{K, h) = l, g x {K, h, hi) =hx+ m dis (K, h) , (12) 

where mdis = rna\ s (K,h) is the magnetization of the disordered state. The latter is thermody- 
namically stable only for T > T c . The surface field hi enters only in the "effective" ordering field 
gi. The so-called extrapolation length A (whose inverse Co = 1/A is conveniently called surface 
enhancement) is positive as it should be if the surface interactions are not enhanced. 

In order to better understand the expression for gi it is helpful to recall some general symmetry 
properties of the Ising Hamiltonian (|l|) which should be respected by the continuum theory. If 
one replaces h and hi by its negative, the mean magnetizations and the local order parameter (^|) 
behave as m; —> —mi, <p n — > — </>„, implying that the ordering surface field should change sign, 
too: 

gi (K,-h,-hi) = -gi(K,h,hi). (13) 

For bulk field h — and arbitrary hi, the Ising antiferromagnet is exactly equivalent to an Ising 
ferromagnet since flipping all Ising spins on one sublattice and changing the sign of K leaves the 
partition function invariant. For the semi-infinite ferromagnet, however, an ordering field gi is 
merely equivalent to a surface magnetic field acting on the spins of the first layer and one easily 
finds go): 

gi(K,0,h 1 ) = hi. (14) 
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The expression (12) fulfills both ( |l3| ) and (|Tj). That the bulk magnetization mdi s of the disordered 
phase comes into play can be understood as follows. For T > T c , the order parameter profile 
should vanish identically if g\{K,h,ht) = 0. Hence the layer magnetization profile m n must 
be a constant, m n = mdi s (K,h). By Eq. Q), the molecular fields acting on surface and bulk 
spins are H + Hi — 4|J| m-i and H — 4 J|(m„_i — rn n+ i), respectively. Thus the flat profile is 
a solution if and only if hi = —m<n s (K,h). Below T c a thermodynamically unstable solution 
m n = 77idj s (i. e. <j) n = 0) still exists if and only if hi = — m^^K, h). In order that this unstable 
solution survives the continuum approximation, one again has to demand that g\(K,h,hi) = 
if hi = -m d i s {K,h). 



4 Wetting phase diagrams 

The physical picture behind the wetting behavior at a (100) surface is the following. Below 
T c , two ordered bulk phases ±</>b coexist. If, e. g., the effective ordering field gi is positive, the 
surface favors the phase <j>b > 0, i. e. A atoms tend to occupy sublattice a planes n = 1, 3, . . . while 
B atoms reside preferentially on the [3 planes n = 2,4, . . .. However, it may occur that deeper 
in the bulk the role of the two sublattices is interchanged and the order parameter assumes the 
value — 4>i, there. Then an antiphase boundary separating two regions of ordered phase appears. 
Antiphase boundaries always occur in real alloys below the ordering temperature. They are the 
analogs of domain walls in an Ising ferromagnet. Now one may ask how the thickness of this 
"adsorbed" layer of bulk phase behaves when the temperature is varied. If the interface stays 
within a finite distance from the surface for T < T w while moving arbitrarily away into the bulk 
for T > T w , a wetting transition takes place at T w — T w (hi, h). The thickness of the layer may 
either grow continuously as T ] T w (second-order wetting transition), or jump from a finite value 
below T w to infinity for T > T w in a first-order transition. 

We have calculated the wetting phase diagram in the space of physical parameters 1/K cx T, 
h, and h\. To this end we employed the full expressions for the bulk and surface Landau free 
energies fi and f s in (J9|) to be found in ||, and used the standard equal-area construction to 
locate the transitions Of course, a description in terms of continuum mean-field theory is 

only sensible above the roughening temperature Tr Jl3|1 , where Tr ~ |T C for h = 1 ||. For 
T < Tr the growth of the wetting layer proceeds via an infinite sequence of layering transitions 
which are outside the scope of the continuum theory. Fig. 2 shows two representative phase 
diagrams at fixed bulk fields h — 1 and h — 1.5. The tricritical point separating continuous 
and first-order transition lines on the right branch of the phase diagram (where gi > 0) is found 
to depend strongly on the bulk field for h larger than « 1.5. If \h\ > 2, the ordered phases are 
energetically unstable at T = and the bulk phase transition ceases to exist. In order to interpret 
Fig. 2 an exact groundstate analysis may be carried out similar to the one in |k| , from which one 
easily finds that complete wetting already occurs at T = if \hi \ > 1. Thus continuum mean- field 
theory is expected to fail if wetting transitions at finite temperatures are predicted for \hi\ > 1. 
In this case one has to resort to the discrete mean-field theory for a more accurate description of 
the low-temperature part of the phase diagram. 
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Figure captions 



Fig. 1. Disordered (A2) and ordered (B2) structures of a bcc binary (AB) alloy. Black and 
white circles denote A and B atoms (or sublattice a and (3 sites), respectively. The hatched 
planes represent the (100) and (110) surface orientations. 

Fig. 2. Wetting phase diagrams at fixed bulk fields h = 1 and h — 1.5, exhibiting continuous 
(full lines) and first-order wetting transitions (dashed lines). One has gi > (gi < 0) to the right 
(left) of the dashed-dotted lines where gi = 0. 
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Figure 1 (Leidl/Drewitz/Diehl) 



(100) (110) 




Figure 2 (Leidl/Drewitz/Diehl) 




